Dynamical study of the singularities of gravity in the presence of 
non-minimally coupled scalar fields 



L.R. Abramo 1 , L. Brenig 2 , E. Gunzig 2,3 , Alberto Saa 2,4 * 
1) Instituto de Fisica, Universidade de Sao Paulo, CP 66318, 05315-970 Sao Paulo, SP, Brazil 

2) RggR, Universite Libre de Bruxelles, CP 231, 1050 Bruxelles, Belgium. 
3) Instituts Internationaux de Chimie et de Physique Solvay, CP 231, 1050 Bruxelles, Belgium. 
4) IMECC - UNICAMP, CP. 6065, 13081-970 Campinas, SP, Brazil. 



Abstract 



We investigate the dynamics of Einstein equations in the vicinity of the 



two recently described types of singularity of anisotropic and homogeneous 



cosmological models described by the action 



S 




with general F(<p) and V(<j>). The dynamical nature of each singularity is 



elucidated, and we show that both are, in general, dynamically unavoidable 



reinforcing the unstable character of previous isotropic and homogeneous cos- 



mological results obtained for the conformal coupling case. 
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I. INTRODUCTION 



We have recently [1] studied the singularities of homogeneous and anisotropic solutions 
of cosmological models described by the action: 

S = J d A x^ {F{<j>)R - d^dy - 2V(cf>)} , (1) 

with general F{<p) and V ((/)). Such singularities appeared in the study of the robustness 
of previously considered [2] homogeneous and isotropic solutions of cosmological models 
governed by (1) with F{4>) = 1 — \4> 2 i corresponding to the so-called conformal coupling, 
and V(4>) = tt0 2 — j0 4 . These homogeneous and isotropic solutions present some novel 
and interesting dynamical behaviors such as: superinflation regimes, a possible avoidance 
of big-bang and big-crunch singularities through classical birth of the universe from empty 
Minkowski space, spontaneous entry into and exit from inflation, and a cosmological history 
suitable in principle for describing quintessence. The appearance of the singularities implies 
that these results are not robust, they are radically changed, even for small disturbances 
in initial conditions and in the model itself. We have shown that the singularities are, 
essentially, of two types. The first one corresponds to the hypersurfaces F((f>) — 0. It is 
not present in the isotropic case, and it implies that all previous homogeneous and isotropic 
solutions passing from the F((f>) > to the F((f>) < region are extremely unstable against 
anisotropic perturbations. The second type of singularity corresponds to Fi((f>) = 0, with 

F 1 (0) = F(0) + ^(F'(0)) 2 , (2) 

and it is present even for the homogeneous and isotropic cases. Although for small deviations 
from the conformal coupling the latter singularities are typically very far from the region of 
interest, in the general case they can alter qualitatively the global dynamics of the model due 
to restrictions that they impose on the phase space. Again, the persistence of some of our 
previously described results, in particular the ones concerning heteroclinic and homoclinic 
solutions, are challenged. 
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Both kinds of singularities have already been described before. To the best of our knowl- 
edge, Starobinski [4] was the first to identify the singularity corresponding to the hypersur- 
faces F((f>) = 0, for the case of conformally coupled anisotropic solutions. Futamase and 
co-workers [5] identified both singularities in the context of chaotic inflation in F{<p) = 1— £</> 2 
theories (See also [6]). The first singularity is always present for £ > and the second one 
for < £ < 1/6. Our conclusions were, however, more general since we treated the case 
of general F((f>) and our results were based on the analysis of true geometrical invariants. 
Our main result is that the system governed by (1) is generically singular on both hyper- 
surfaces F((f>) = and Fi(<f>) = 0. Here, generically means that it is possible to construct 
non-singular models if one fine-tunes F(4>) and V((f>), as we have shown in [1]. The physical 
relevance of such a fine-tuned model is still unclear. 

As was shown in [1], one can advance that there are some geometrically special regions in 
the phase space of the model in question by a very simple analysis of the equations derived 
from the action (1). They are the Klein-Gordon equation 

□0 - V'(<f>) + l -F'{<t>)R = 0, (3) 

and the Einstein equations 

F{<f>)G ab = {l + F"{4>))d a 4>d b 4> 

- l -g ah [(1 + 2F"(0))9 C 09 C + 2V {<!>)] - F'{4>) {g ab U4> - V a 0V 5 0) . (4) 

We consider now the simplest anisotropic homogeneous cosmological model, the Bianchi 
type I, whose spatially flat metric is given by 

ds 2 = -dt 2 + a 2 (t)dx 2 + b 2 (t)dy 2 + c 2 (t)dz 2 . (5) 

The dynamically relevant quantities here are 

H, = -, H 2 = b -, andtf 3 = -. (6) 
a b c 

For such a metric and a homogeneous scalar field <fi = (f>(t), after using the Klein-Gordon 

Eq. (3), Eq. (4) can be written as 
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F(<f>)G 00 = \ft + V{4>) - F\4>) (H, + H 2 + H 3 ) 0, (7) 
^)Gn = 1 + 2 f (0) 2 " ^(0) - F'(0) ^ + V(4>) - ^i?) , (8) 
^F{<t>)G 22 = 1 + 2 f (0) 2 - Vtf) - F\<t>) [h^ + V'W - ^-R^j , (9) 

1f(0)g 3 3 = 1 + 2 f (0) 2 - v(0) - f\j>) [h^ + \/'(0) - • (10) 

It is quite simple to show that Eqs. (8)-(10) are not compatible, in general, on the hy- 
persurface F(<f) = 0. Subtracting (9) and (10) from (8) we have, on such hypersurface, 
respectively, 

F^Xifi - H 2 )j> = 0, and F'^H, - H 3 )<j> = 0. (11) 

Hence, they cannot be fulfilled in general for anisotropic metrics. As it was shown, this 
indeed corresponds to a geometrical singularity which cannot be prevented in general by 
requiring that F'{<p) = or = on the hypersurface. 

As to the second singularity we have, after taking the trace of the Einstein equations, 
that: 

R = 0) = (4^(0) + 3\/'(0)F'(0) - (1 + F"{4>))4> 2 ) ■ (12) 

Inserting Eq. (12) in the Klein-Gordon Eq. (3), one can see that it contains terms which are 
singular for Fi((f>) = 0. Again, as we will see, this corresponds to an unmovable geometrical 
singularity, and it cannot be eliminated, in general, by demanding that F'(<f>) = on the 
hypersurface Fi(<f>) = 0. On both the hypersurfaces F((f>) = and Fi((f>) = the Cauchy 
problem is ill-posed, since one cannot choose general initial conditions. 

The hypersurfaces F(<f) = and F\{4>) — also prevent the global definition of an 
Einstein frame for the action (1), defined by the transformations 

9ab = F((f))g ab , (13) 

(dJ\ 2 = FM_ 



It is well known that in the Einstein frame the Cauchy problem is well posed. Again, 
the impossibility of defining a global Einstein frame shed some doubts about the general 
Cauchy problem. Moreover, the standard perturbation theory for helicity-2 and helicity-0 
excitations, derived directly from Eqs. (13)-(14), fails on both hypersurfaces [7]. 

The question to be addressed in the following sections is the dynamical behavior of the 
Eqs. (3) and (7)-(10) in the vicinity of the two hypersurfaces corresponding to F((j)) = 
and = 0. As we will see, both hypersurfaces are dynamically unavoidable, meaning 

that they have an attractive neighborhood, excluding definitively the possibility that these 
singularities are hidden by some dynamical barrier that would prevent the solutions to reach 
them. Whenever a solution enter in the attractive neighborhood, it will unavoidably reach 
the singular hypersurface. 

II. NON-CONSERVATIVE SYSTEMS AND THE DIVERGENCE THEOREM 

For Hamiltonian systems, as a consequence of Liouville theorem, phase space volumes are 
preserved under the system time evolution. That means that if one chooses an initial closed 
hypersurface So in the phase space and let each point of S evolve in time according to the 
system equations, the closed hypersurface S will evolve to another closed hypersurface S t 
at some latter time t, and the volumes V of the region enclosed by S and S t are exactly the 
same, V(0) = V(t). This is a characteristic of conservative systems. The system given by 
(3) and (7)-(10) is not conservative. By choosing the set of coordinates (0, ip,p, q, r) (with 
ip — (j), p — H 1 + H 2 + H 3 , q = Hi — H 2) and r = Hi — H 3 ), for the phase space V, Eq. (3) 
and (7)-(10) can be cast in the form 



(<j),ip,p,q,f) = W((f),ijj,p,q,r). 



(15) 



For the metric (5), we have the following identities 



Coo — HiH 2 + H 2 H 3 + HiH 3 , 
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G 22 = b 2 [k 2 + H^H, + H 2 + H 3 ) - ^R) , (16) 
G33 = c 2 (h 3 + H 3 (H X + H 2 + H 3 ) - ^R) , 
R = 2(h 1 +H 2 + H 3 + Hi + Hi + Hi + H X H 2 + H 2 H 3 + ifxifg) • 

— * 

Using them, the components of W can be explicitly computed from Eqs. (7)-(10), 
= -jnl> - V\<f>) + ^F\<f>)R^), 



W p 



(F(0) + 2F'(<f>) 2 )p 2 + |(1 + 2F"(0))0 2 - W{4>) - 3F'(0)V'(0) 
+ (F(0) + F'(0) 2 )(g 2 + r 2 - gr)] /(2F 1 (c/ ) )) 

Wr= -(p + ™*y. (17) 

The divergence theorem assures us that the volume V of a closed hypersurface of P 
evolves in time as: 

V(t)= I (divW')dvol, (18) 

J St 

where the integral is performed in the region enclosed by St. The divergence of the vector 

— * 

field W determines, therefore, how fast a volume of a closed hypersurface of V is expanded 
(diviy > 0) or contracted (divjy < 0). For conservative systems, one has divVV = 0. A 
straightforward calculation here give us 

div W = -p - 2 f p + ^J^-iA 



(F(0) + 2F'(0) 2 ) p + (f'(0) (1 + F"{4>)) - \F'{4>)) i> 



(19) 



It is clear from (19) that our system suffers violent contractions and/or expansions in the 
neighborhood of the hypersurfaces F{<p) = and -F\(0) = 0. Let us consider each of them 
separately since they lead to different kinds of singularity. 
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III. PHASE SPACE CONTRACTION AND EXPANSION NEAR F(<f>) = 



Homogeneous and isotropic solutions, for whose q = r = 0, are known to be perfectly 
regular on F{(f>) = [2], in contrast to (19) that presents unequivocally a divergence on 
this hypersurface. A closer analysis of the vector field W (17) reveals that the divergent 
contractions and expansions near F((fi) = are associated to the directions q and r, namely 
the quantities that measure the anisotropy of the solution. The other directions of the 
flux defined by W are regular on the hypersurface F(<f>). That means that if a solution is 
perpendicular to the divergent directions q and r, solutions for which q = r = 0, i.e. isotropic 
solutions, it will evolve without suffering any violent contraction or expansion in the other 
directions of V . Since W q and W r are proportional to q and r, respectively, an initially 
isotropic solution g(0) = r(0) = remains isotropic for all latter t, q(t) = r(t) = 0. We can 
say that isotropic solutions are orthogonal to the divergent fluxes. However, any amount, 
no matters how small, of anisotropy (non vanishing q or r) will break the orthogonality and 
the solution will prove the divergent directions, being strong contracted or expanded. This 
is the dynamical origin of the instabilities of anisotropic solutions near F(4>) = 0. 

Let us suppose now that F'((p) ^ on the hypersurface F((p) = (if F'(<f>) vanishes on 
the hypersurface F'(<f)) = 0, then by Eq. (2) both hypersurfaces F(<p) = and Fi(<f>) = 
coincide). The corresponding pole on O (F(0 O ) = 0) in the volume integral (18) will have 
as numerator the factor — 2F'(0 o )'*/s implying that the flux defined by (17) passes from a 
catastrophic contraction to a catastrophic expansion as one passes by 4>q. Since = ip (see 
Eq. (17)), any solution approaching the hypersurface F(<p) = with ip ^ (we excluded 
from the analysis the possibility of having fixed points on F((f>) = 0, for these cases, of 
course, it makes no sense to talk about "crossing" F(4>) = 0) will cross it and, hence, prove 
the divergent phases of contraction and expansion. 

In the expanding "side" of the hypersurface F((p) = 0, q and r diverges as — > (fio, and 
the system will be unavoidably driven toward a spacetime singularity [1], as we can conclude 
by considering, for instance, the Kretschman invariant / = R abcd R abcd , which for the metric 
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(5) is given by 

/ = 4 + H{f + (H 2 + iJ|) 2 + (H 3 + iJ 3 2 ) 2 + tf 2 tf 2 2 + HlHl + tf 2 tf 3 2 ) . (20) 

The invariant / is the sum of non negative terms. Moreover, any divergence of the variables 
Hi, H 2 , H 3 , or of their time derivatives, would suppose a divergence in J, characterizing a 
real geometrical singularity. Since the relation between the variables p, q, r, and Hi, H 2 , H 3 
is linear, any divergence in p, q, r or of their time derivative, will suppose a divergence in 1. 

IV. PHASE SPACE CONTRACTION AND EXPANSION NEAR Fi{4>) = 

The hypersurface Fi(<f>) = 0, as the previously one considered in the last section, also 
separates regions of catastrophic contraction and expansion in the phase space V . However, 
in contrast to the previous one, this hypersurface leads to singularities even for homogeneous 
and isotropic solutions. The divergent directions now are p and ip, and there are no solutions 
of (15) orthogonal to them, since, in contrast to the q and r directions, for which W q and 
W r are respectively proportional do q and r, and W p do not vanish for ip = and 
p = 0. In this case, no solution can escape from crossing = 0. In the expanding side 

of Fi(<f>) = 0, p — > oo as <j) — > 0i (Fi(0i) = 0), implying the divergence of the invariant (20). 

V. FINAL REMARKS 

The singularities described in the precedent section imply that the model presented in 
[2,3] is not robust, since its main conclusions were a consequence of very especial initial 
conditions, i.e. they are valid only for solutions orthogonal (namely the isotropic ones) to 
the divergently expanding directions q and r. For instance, all homogeneous and isotropic 
solutions crossing the F((f>) = hypersurface are extremely unstable against anisotropic 
perturbations. Any deviation from perfect isotropy (expressed by nonvanishing q and r 
variables) for these solutions, however small, will lead catastrophically to a geometrical sin- 
gularity. Many of the novel dynamical behaviors presented in [2,3] depend on these solutions. 
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This is the case, for instance, of some solutions exhibiting superinflation regimes. The het- 
eroclinic and homoclinic solutions identified in [2,3] can cross the F{<p) = hypersurface 
and, hence, they also suffer the same instability against anisotropic perturbations. The 
homoclinic solutions were considered as candidates to describe a non-singular cosmological 
history, with the big-bang singularity being avoided through a classical birth of the universe 
from empty Minkowski space. Apart from F((f>) = singularities, these solutions are also 
affected by the singularities of the type Fi(<f>) = 0. Suppose that the conformal coupling is 
disturbed by a very small negative term: F{<p) — 1 — (| — e)0 2 . The Fi(<j>) = singularities 
will be near the = ±l/y/e hypersurfaces. Although they are located far from the F(</>) = 
regions, they alter the global structure of the phase-space. In this case, they restrict the 
existence of homoclinics, rendering a non-singular cosmological history more improbable. 

The singularities do not affect the conclusions obtained by considering solutions inside 
the F(4>) > region. The asymptotic solutions presented in [3], for instance, are still valid. 
The conclusion that for large t the dynamics of any solution (inside F(<p) > 0) tends to an 
infinite diluted matter dominated universe remains valid. Moreover, for small anisotropic 
deviations (q and r small in comparison with p), the solutions inside F(<f>) > 0, for large t, 
approach exponentially isotropic matter-dominated universe. 

A singularity-free model can be constructed by demanding a well behaved divW on 
both hypersurfaces. This can be achieved [1] by requiring F(0 O ) = F'(0 O ) = 0, and by 
choosing a V{4>) that goes to at a proper rate when — > <p . Moreover Fi(0) must 
have no other zeros than the ones of F(<f)). Models for which F(4>) = (<p 2n and V(<f>) = 
a( p2(2n-i) _|_ hjgh or d er terms, for instance, fulfill these requirements. However, such a highly 
fine-tuned class of model is of no physical interest here, since it does not contain F{4>) > 
and F(4>) < regions and consequently has no solution for which the effective gravitational 
constant G c fj changes it sign along the cosmological history. The stability of such solutions 
were the starting point of the analyses of the pioneering work [4] and of the present one as 
well. 
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